Observational tests of inflation with a field derivative coupling to gravity 
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A field kinetic coupling with the Einstein tensor leads to a gravitationally enhanced friction during 
inflation, by which even steep potentials with theoretically natural model parameters can drive 
cosmic acceleration. In the presence of this non-minimal derivative coupling we place observational 
constraints on a number of representative inflationary models such as chaotic inflation, inflation 
with exponential potentials, natural inflation, and hybrid inflation. We show that most of the 
models can be made compatible with the current observational data mainly due to the suppressed 
tensor-to-scalar ratio. 



I. INTRODUCTION 

Inflation has been the backbone of the high-energy cos- 
mology over the past 3 decades [l|. The most simple 
source for inflation is a minimally coupled scalar field 6 
("inflaton") with a slowly varying potential V(4>) @, HJ 
The spectra of density perturbations generated from the 
quantum fluctuations of inflaton are consistent with the 
temperature anisotropics observed in the Cosmic Mi- 
crowave Background (CMB) j|. 

From the amplitude of the observed CMB anisotropics 
0, 0] the typical mass scale of inflation is known to be 
around m ~ 10 14 GeV Q. This is much larger than 
the electroweak scale (~ 10 2 GeV), which suggests the 
requirement of new physics beyond the Standard Model 
of particle physics [8j. In other words, for the potential 
V((j>) = (A/4)(</> 2 - v 2 ) 2 with v ~ 10 2 GeV, the coupling 
A is constrained to be A ~ 10~ 13 from the CMB normal- 
ization Q, but this is much smaller than the coupling 
constant A ~ 0.1 of the Higgs boson 

There have been attempts to accommodate the Higgs 
field for inflation. One is to use a non-minimal field 
coupling £,R4> 2 1 2 with the Ricci scalar R [To[ (see also 
Refs. [Hi). If £ S> 1 the self coupling A can be as large 
as A w 10" 10 £ 2 from the CMB normalization [ij|. Al- 
though this scenario is attractive, it is plagued by the 
unitary-violation problem associated with graviton ex- 
change in 2 —> 2 scalar scattering around the energy scale 
A c « M p i/£ (where M pl = 2.44 x 10 18 GeV is the reduced 
Planck mass) [li|. Since A c is around the energy scale 
of inflation, some strong coupling effect can give rise to 
additional corrections to the inflaton potential. 

Another attempt is to employ a field deriva- 
tive coupling with the Einstein tensor G^ v ', i.e. 
G^ v d a 4)d v <t) / (2M 2 \ where M is a constant having a di- 
mension of mass [l4| ( see also Ref. [HI for the original 
work). In the regime where the Hubble parameter H 
is larger than M the field evolves more slowly relative 
to the case of standard inflation due to a gravitationally 
enhanced friction. Hence it is possible to reconcile steep 
potentials such as V{<j>) = A0 4 /4 (A ~ 0.1) with the CMB 
observations. 

In Refs. [l3, Il6l - fl9j it was shown that, for a 
slow-rolling scalar field satisfying the condition e = 



{d(f>) 2 1 \M 2 M 2 X ) < 1, the strong coupling scale A c of the 
derivative coupling theory is around M p \ in a homoge- 
neous and isotropic cosmological background. Provided 
that H and M are below the Planck scale, the theory 
is in a weak coupling regime with suppressed quantum 
corrections. 

The property of the high cut-off scale A c around M p \ is 
associated with the fact that whenever the non-minimal 
derivative coupling to gravity dominates over the canoni- 
cal kinetic term the theory possesses an asymptotic local 
shift symmetry for e < 1 [19]. This symmetry is re- 
lated to the Galilean symmetry <f> — > <fi + c + c^x 11 in 
Minkowski space-time [20] , but the difference is that the 
coordinate x^ in the derivative coupling theory on curved 
backgrounds is linked to the covariantly constant Killing 
vectors [l!|. In the presence of a slowly varying inflaton 
potential such a local symmetry is only softly broken, so 
that the potential can be protected against quantum cor- 
rections during inflation. The field self-interaction of the 
form (d(f) 2 D(j) pol l2lj. which satisfies the Galilean sym- 
metry in the limit of Minkowski space-time, also leads to 
the slow evolution of <fr along the inflaton potential [22[ 
(see also Refs. H3]). 

A nice feature of the non-minimal derivative coupling 
with the Einstein tensor is that the mechanism of the 
gravitationally enhanced friction works for general steep 
potentials. For instance, let us consider the potential 
of natural inflation, V((j>) = A 4 [l + cos(0//)], where / 
characterizes the scale of the breaking of a global shift 
symmetry [24| ■ In order for this potential to be con- 
sistent with the CMB observations, wc require that f i s 
larger than 3.5M p i in conventional slow- roll inflation [251 ] . 
Then the global symmetry is broken above the quantum 
gravity scale, in which case quantum field theory may be 
invalid [26j . In the presence of the non- minimal deriva- 
tive coupling to gravity, however, the scale / can be much 
smaller than M p \ because of the gravitationally enhanced 
friction [11 [3. 

Another example is the exponential potential V(<fi) = 
y Qe P<f>/M p i ^ wnose dominance leads to the power-law ex- 
pansion of the Universe (with the scale factor a oc t 2 ^ , 
where t is cosmic time) [271 ] . In higher-dimensional grav- 
itational theories, exponential potentials often arise as 
the curvature of internal spaces related with the geome- 
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try of extra dimensions [28[. In such cases the constant 
/? is usually larger than the order of unity, so that it 
is difficult to realize sufficient amount of inflation. As 
we will see later, this problem can be circumvented by 
taking into account the non-minimal derivative coupling. 
Moreover, unlike the standard case, inflation comes to 
end with gravitational particle production. 

In order to test the viability of inflationary models with 
the non-minimal derivative coupling it is important to 
estimate the power spectra of density perturbations rel- 
evant to the CMB anisotropics. In Refs. [Ill [l8[ the 
authors computed the inflationary observables such as 
the scalar spectral index n s , the tensor-to-scalar ratio 
r, and the nonlinear parameter /nl"' 1 of the eq uilateral 



scalar non-Gaussianities (see also Refs. |29|-[31|). Since 
the scalar propagation speed is close to the speed of light 
during inflation, the scalar non-Gaussianities are sup- 
pressed to be small (I/nl" 1 '! ^ !)■ Hence n s and r are 
the two main observables to distinguish between different 
inflaton potentials. 

In this paper we shall place observational constraints 
on a number of representative inflationary models in the 
presence of the field derivative coupling with the Ein- 
stein tensor. We use the bounds derived from the joint 
data analysis of WMAP7 @, Baryon Acoustic Oscilla- 
tions (BAO) [32|, and the Hubble constant measurement 
(HST) [33[. Note that some constraints on Higgs inflation 
and natural inflation have been discussed in Refs. [l6l.H8j 
without the CMB likelihood analysis. In Ref. [34| the au- 
thor carried out the cosmological Monte-Carlo simulation 
to test Higgs inflation with the field derivative coupling to 
gravity. Our analysis based on the recent observational 
data is general enough to cover a wide variety of mod- 
els such as chaotic inflation Q , inflation with exponential 
potentials p7| , natural inflation [24[ , and hybrid inflation 
[3a ] . We show that the gravitationally enhanced friction 
mechanism can make most of the models compatible with 
the current observations. 



II. BACKGROUND DYNAMICS 



We start with the following 4-dimcnsional action 



S = 



where 



Kl 1 



,<t>-v{4>) 



M 2 



(1) 



(2) 



Here g is a determinant of the space-time metric g^ v , R 
is the Ricci scalar, G^ u is the Einstein tensor, M is a 
constant having a dimension of mass, and V((f>) is the 
potential of a scalar field (p. 

The action ([]} belongs to a class of the most gen- 
eral scalar-tensor theories having second-order equa- 
tions of motion (which is required to avoid the Os- 
trogradski instability) [2|| l3q - l38j . The Lagrangian 



in such general Horndcski's theories is the sum of 
the terms C 2 = K{<j>,X), £3 = -G 3 (</>, A)D0, 
£4 = Gi((j),X)R + G41 x [field derivative terms] , 
and C 5 = G 5 (0,A)G^(V Al V^) - (G 6rX /6) x 
[field derivative terms], where K, Gi (i = 3,4,5) are 
functions of <fi and X = — g^ v d^4>d u (j)/2 1 and Gi.x = 
dGi/dX [2^, H3]. The conditions for the avoidance 
of ghosts and Laplacian instabilities were derived in 
Refs. [2^, [39]]. These conditions can be used to restrict 
the functional forms of K, Gi (i = 3, 4, 5) to construct 
theoretically consistent models of inflation. 

The non-minimal derivative coupling in Eq. ([T]) is re- 
covered in the Horndeski's Lagrangian by choosing the 
function G5 = —(f)/ (2M 2 ) after integration by parts. The 
sign in front of the term G^ v /M 2 in Eq. @ is chosen 
to avoid the appearance of ghosts in the scalar sector 

HE!. 

In Ref. it was found that in a manifold hav- 
ing integrable (covariantly constant) Killing vectors £ a 
the field Lagrangian — A A " y 9 M 09 I/ (/)/2 in Eq. |T]) is in- 
variant under the (curved-space) Galilean transformation 
4>(x) — > <f)(x) + c + c a £ a , where c, c a , xo are constants 
and a; is a space-time coordinate. The existence of the 
Galilean symmetry has an advantage that the theory can 
be quantum mechanically under control [4(1 ] . 

Imposing the above Galilean symmetry in the curved 
background with integrable Killing vectors, Gcrmani et 
al. showed that the second-order Lagrangians are re- 
stricted to take the forms ~A fJ, ^d^4>d 1 y(j)/2 or XOcj) (plus 
a field derivative coupling with the double dual Riemann 
tensor). In the small derivative regime in which the 
condition (d(j)) 2 / (M 2 M 2 X ) < 1 is satisfied (e.g., during 
inflation), an approximate infinitesimal shift symmetry 
4> — > <f) + f{x) (where f{x) is an arbitrary function of 
space-time coordinates x) emerges for the Lagrangian 
G^d fl cf)d u (t)/(2M 2 ), provided that the metric is shifted 
appropriately [H^ . The existence of such a gauge sym- 
metry can allow the theory ([T]) to be protected against 
quantum corrections even up to the Planck scale. Note 
that the term XOcj) does not possess such a general gauge 
shift symmetry. 

The scale of unitarity violation for the theory (JTJ) was 
estimated in Ref. [l4| in the context of Higgs inflation. 
In Standard Model we can consider the <j)(f) — > (jxj> scat- 
tering via graviton exchange, where 4> is one of the real 
scalar degrees of freedom for the Higgs doublet. We 
expand the metric in the form g M „ = g$ + h^ v /M p \, 

where gffl = (—l,a 2 (t),a 2 (t),a 2 (t)) is the metric on 
the flat Friedmann-Lemaitre- Robertson- Walker (FLRW) 
background (a(t) is the scale factor with cosmic time t). 
We are interested in the high-friction regime in which the 
Hubble parameter H = a/ a (a dot represents a derivative 
with respect to t) is much larger than M. In this regime 
the field <f> is expanded as <f) = 0o + M\/ (>/3H), where x 
is a canonically normalized field perturbation. The first 
non-renormalizablc operator associated with the interac- 



3 



tion between gravitons and scalars is given by \\A 



2H 2 M pl 



(3) 



A power counting analysis gives the unitary bound A ~ 
{2H 2 M pi ) 1 / 3 . For the suppression of higher dimensional 
operators we require the condition R < A 2 . On using the 
relation R ~ 12H 2 this condition translates into H < 
5 x 10 _2 M p i, which is satisfied during inflation. 

The discussion of the unitary bound given above can 
be applied to the multi-field inflationary models in which 
one of the fields is not necessarily responsible for the 
cosmic acceleration. In single field models where only 
one field (f> leads to inflation the unitary bound can be 
as close as the Planck scale M p \ in the regime where the 
condition {d(j)f / (M 2 M^) < 1 is satisfied [3 EH- In 
this case the slow-roll evolution of the field <f> suppresses 
the interaction ([3]) below the Planck scale. 

In the following let us study the background dynam- 
ics for the theory described by the action ([IJ. In the 
flat FLRW background the equations of motion follow- 
ing from the action ((T|) are 



where 



A = l 



•EL 

'm 2 



(13) 



We define the following slow-roll parameters 



H 
IP 



Sx 



2H 2 M 2 l ' 



H4> 
Sd = 



AM 2 M 2 



(14) 



pi 



For the validity of the slow-roll approximation we require 
that {e, \S,p\,Sx,Sd} -C 1. Taking the time-derivative of 
Eq. (fTTj) and using Eq. (TT2"j) . we have 



'A' 



where 



M 2 pl 



Ya 

V 



(15) 



(16) 



SM^H 2 = 



l__d 



Hi 

V 2 







a 3 







+ Va = : 



(4) 
(5) 



where = dV/d<f>. In order to solve the dynamical 
equations numerically, it is convenient to introduce the 
following dimcnsionless variables 



M p i 



y 



MM p i 



H 

M 



(6) 



Differentiating Eq. Q with respect to t and using Eq. ([5]) 
to eliminate H, we obtain the second-order equation for 
the field 6. It then follows that 



dx 

dy 
<1t 



where 
t = Alt. 



3yz(2 - 3y 2 )(l + 3^ 2 ) + (2 - y 2 )V^{x) 
2(l + 30 2 ) + y 2 (9z 2 - 1) 



! y 2 + 2V{x) 
6 - 9y 2 



V, 



M 2 M, 



V 



V 



pi 



M 2 M 2 vl 



(7) 

(8) 

(9) 
(10) 



We are interested in slow-roll inflation in which cos- 
mic acceleration is mainly driven by the potential energy 
V{4>). In this case Eqs. dU and ([5]) reduce to 



UI^H 2 



V{4>), 
3HAc/) + V (j> ~ 0. 



(11) 
(12) 



This shows that e <C e v for A ^> 1 and hence the evolu- 
tion of the field (f> slows down relative to that in standard 
slow-roll inflation. 

The field value <pf at the end of inflation is known by 
solving e((f>f) = 1, i.e. 



£ v(^/) 



M 2 M 2 X 



1 . 



(17) 



The number of c-foldings from the time t during inflation 
to the time tf at the end of inflation is defined by N = 
Jl s H(t) dt. On using Eqs. CP and @2J), it follows that 



N 



Ml 



V 



4> f 



V 



M 2 M 2 J V$ 



(18) 



If A 3> 1 (i.e. H 2 3> M 2 ) during inflation, we can neglect 
the first term inside the bracket of Eq. (|18|) relative to the 
second one. This is not the case for inflation in which the 
transition from the regime H > M to the regime H < AI 
occurs prior to the onset of reheating. 

As an example, let us consider chaotic inflation Q with 
the potential 



v{4>) 



X 



(19) 



where A and n are constants. In this case Eqs. (fTS"]) and 
(fI7]l read 



N 



2x) 



x 
2n 



2 r 



2o 



( 1 + ^?) 



n(n + 2) 

,2 



2n 



2a 



n(n + 2) 



,(20) 
(21) 



4 



s.o 




condition is violated. In this regime some quantum cor- 
rections may come into play to the action ([TJ . As long as 
such corrections are unimportant in the field equations 
(U]) and (JSJ) , we find that the inflaton oscillation is not 
disturbed during the transient period (see Fig.Q}. In the 
case (ii) we have zt = 0.3 at Xf = 1.3 and hence the 
system enters the regime H < M during inflation. In 
this case the qualitative behavior for the oscillation of 
inflaton at reheating is not much different from that in 
standard inflation. 

We confirmed that the difference between the numer- 
ical and analytic values of the number of e-foldings ac- 
quired during inflation is usually less than a few percent. 
This shows that the slow-roll approximation employed in 
Eqs. (fTTT) and (|T2"j) can be trustable. 
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Figure 1: Evolution of the field cj> versus r = Mt for (i) a — 
100 and (ii) a = 0.25. The initial conditions at iV = 60 are 
chosen by using the slow-roll equations {TTJ and <|12p . i.e., (i) 
x = 1.757, y = -0.526, and (ii) x = 7.397, y = -0.521. 



where 



AM 



n-2 
Pi 



A/ 2 



For the quadratic potential V{4>) = my/2 (i.e. A = m 
and n = 2) one has x 2 = + 4a — l)/a and 

x 2 = ^y / 2 + 2Vl + 4a + 4a(4AT + l) -2^ /a, (23) 

where a = m 2 /M 2 . In the General Relativistic (GR) 
limit (a — > 0) this gives x 2 — > 4N + 2. In the high- 
friction limit (a — > oo) one has x 2 — > 2^/ (4N + l)/a, 
which means that the field value is smaller than that in 
standard chaotic inflation. 

In order to confirm the accuracy of the slow-roll ap- 
proximation we solve the full equations of motion ©-(O 
numerically. In Fig. Q] the evolution of the field (j> is plot- 
ted for the potential V(<f>) = m 2 <f> 2 /2, i.e. V$ = ax and 
V = ax 2 /2 in Eqs. (0 and ©. We choose the initial 
conditions of x and y at N = 60 by using the values 
derived under the slow-roll approximation. 

The numerical simulations labeled as (i) and (ii) in 
Fig. [1] correspond to the parameters a = 100 and a = 
0.25, respectively. In the case (i) the numerical value 
of z = H/M at the end of inflation (xf = 0.44) is 
Zf = 2.2, which means that the solution is in the high- 
friction regime during inflation. In this case H drops 
below M at the reheating stage. After inflation there is 
a transient period with H > M in which the slow-roll 



III. THE SPECTRA OF DENSITY 
PERTURBATIONS 

The spectra of scalar and tensor perturbations gener- 
ated in the theories given by the action (|T|) were derived 
in Refs. 0, [TH, H^, [3l[. Here we briefly review their 
formulas in order to apply them to concrete inflationary 
models in Sec. IIVI 

The perturbed metric about the flat FLRW back- 
ground is given by [4l| 

ds 2 = -(l + 2A)dt 2 + 2diBdtdx l 

+a 2 (t) [(1 + 2K)8 tJ + hij] dx l dx j , (24) 

where A, B, and 1Z are scalar metric perturbations, and 
hij are tensor perturbations which are transverse and 
traceless. The spatial part of a gauge-transformation 
vector is fixed by gauging away a perturbation E ap- 
pearing as a form Ey in the metric (f2"4")) . We decompose 
the inflaton field into the background and inhomogeneous 
parts, as </> = <j>a{t) + Sc/>(t, x). In the following we choose 
the uniform-field gauge characterized by 5(f) = 0, which 
fixes the time-component of the vector 

Expanding the action (JTJ) up to second order in per- 
turbations and using the Hamiltonian and momentum 
constraints, we obtain the second-order action for scalar 
perturbations [1^, |3l| 



dt d 3 xa 3 Q s 



n 2 



where 



Q s = 



Wi(Aw\W3 + 9W 2 .) 
3^1 



(25) 



(26) 



2 3(2lV 2 W2H — W%W4 + A1V1W1W2 — 2w 2 W2) 



and 



Wl = M 2 , (1-25 D ), w 2 = 2HM 2 l (1-65 D ). 
w 3 = -3fT 2 M£(3 -Sx- 368 D ) , 



W4 = Af5(l + 25 D ) . 



(28) 
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In order to avoid the appearance of scalar ghosts and 
Laplacian instabilities we require that Q s > and c 2 > 
0. Picking up the dominant contributions to Q s and c 2 
under the slow-roll approximation, we obtain 



M^(Sx + 65 D ) 

25 D {3S X + M6 D 



pl A ' 



1 



2<5 ) 



65 



D 



(29) 
(30) 



which mean that c 2 = 1 — 0(e). The power spectrum 
of the curvature perturbation 1Z, which is evaluated at 
c s k = aH (k is a comoving wavenumber), is given by 



V„ = 



H 2 



V 3 



V 



M 2 M 2 p 



(31) 



where in the last approximate equality we used Eqs. ([lip . 
(PI). OH), (HU), and O- The scalar spectral index is 



n H — 1 = 



din A: 
1 

~A 



c s k—aH 



2e v 4 



A 



where e v is dchned in Eq. (|T6|) . and 



Vv 



V 



(32) 



(33) 



- 1 



In the high-friction limit (A S> 1) one has n s 
-(8e v -2r] v )/A with A ~ V/(M 2 M^), whereas n s -l ~ 
-6e v + 2rj v in the GR limit (A ~ 1). 

The intrinsic tensor perturbation /i^ can be decom- 
posed into two independent polarization modes, i.e. 
hij = h + e~lj + h x e*. In Fourier space we normalize the 
two modes, as ef -(fe)ef-(- fe)* = 2 (where p = +, x) and 

efAk)efA— k)* = 0. Then the second-order action for 

m i — h 

tensor perturbations can be written as [29J, ^1[ 



which is evaluated at c t k 
index is 

dlnVt 



aH. The tensor spectral 



-2c 



c t k—aH 



16e v - 



d\n k 

The tensor-to-scalar ratio is 

^ s k~aH 

from which we obtain the consistency relation 
r ~ — 8nt . 



.4 



16e. 



(38) 



(39) 



(40) 

This relation is the same as that in conventional inflation 
at leading order in slow-roll. 

The running spectral indices a s = dn s /dhi k\ Cs k =a H 
and at = dnt/d\nk\ Ct k= a H are second order in slow-roll 
parameters. They are set to be in the CMB likelihood 
analysis. The consistency relation (|40| reduces the infla- 
tionary observables to three, i.e., n s , r, and V s . These 
observables are varied in the likelihood analysis with the 
pivot wavenumber fco = 0.002 Mpc -1 , by assuming the 
flat A-cold-dark-matter model. 

Since c 2 = 1 — 0(e), the general formula for the equilat- 
eral non-Gaussianities of scalar perturbations (3l| shows 
that the nonlinear parameter /nl" 1 is °f the order of e 
fl8| . Hence the scalar non-Gaussianities do not provide 
additional constraints to those derived from the linear 
perturbations. 



IV. OBSERVATIONAL CONSTRAINTS 

In the presence of the field derivative coupling to the 
Einstein tensor we place observational constraints on a 
number of models such as (i) chaotic inflation, (ii) infla- 
tion with exponential potentials, (hi) natural inflation, 
and (iv) hybrid inflation. Our analysis covers most of 
the representative inflaton potentials proposed in litera- 
ture. 



S, 



(2) -E 
p 



dt d 3 xa 3 ( 



where 



Qt 

■2 



wi/4 = 



a 1 



28 D )/A, 
4<^ + 0(e 2 ) 



(34) 



(35) 
(36) 



This shows that, unlike the scalar propagation speed 
squared, c 2 is slightly larger than 1 during inflation. Since 
Lorentz invariance is explicitly broken on the FLRW 
background, the superluminal mode does not necessarily 
imply a violation of causality. The tensor power spec- 
trum is given by 



H 2 



2V 



2n 2 Q t ct 



3^M p 4 



(37) 



A. Chaotic inflation 

We start with chaotic inflation characterized by the 
potential 

Vffl = V , (41) 
n 

where, for n = 2, we use the notation A = m 2 as in the 
previous section. Under the slow-roll approximation the 
dimensionless field x = <f>/M p i is related with the number 
of e-foldings TV as Eq. From Eqs. and (JSgJl the 

scalar spectral index and the tensor-to-scalar ratio are 
given, respectively, by 

- 2 [n(n + 2) + 2(n+l)aa; n ] 



= 1 - 



x 2 (n + ax n ) 2 



8n s 



x 2 (n 



ax 



r 



(42) 



(43) 
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Figure 2: Observational constraints on chaotic inflation with 
the quadratic potential V((j>) = m 2 (f> 2 /2 in the (n B ,r) plane 
with the three different numbers of e-foldings (N = 50, 60, 70). 
We evaluate the theoretical values of n s and r in the range 
10~ 8 < a = m /M 2 < 10 8 . The thin solid curves correspond 
to the lcr (inside) and 2a (outside) observational contours 
constrained by the joint data analysis of WMAP7, BAO, and 
HST. For larger values of a = m 2 /M 2 the tensor-to-scalar ra- 
tio r gets smaller, whereas the scalar spectral index increases. 



Figure 3: Observational constraints on chaotic inflation with 
the quartic potential V(<t>) = X(f> 4 /A in the range 10 -8 < a = 
\{M p i/M) 2 < 10 8 with the three different values of N. The 
lcr and 2cr observational contours are the same as those in 
Fig. [2] While the standard case (a = 0) is outside the 2a 
bound, the field derivative coupling to gravity can make the 
quartic potential compatible with observations. 



2 /M 2 . 
In the limit that a 



where a = XM™ { 



one has x = v / 2n(4N + n)/2 
from Eqs. (f2(J| and (j!?Tj) . in which case n s and r are 



1 



2(n + 2) 

~ AN + n ' 
16n 



AN- 



n 



(44) 
(45) 



These values correspond to those for standard chaotic 
inflation. If N = 60, then n s = 0.967, r = 0.132 for 
n = 2 and n s — 0.951, r = 0.262 for n — A. In another 
limit a — > oo one has x n+2 = [2N(n + 2) + n]n 2 /(2a) and 



7l s = 1 - 



4(71+1) 

2(n + 2)N + n 
16n 



2(n + 2)N + n 



(46) 
(47) 



If N = 60, then n s = 0.975, r = 0.066 for n = 2 and 
n s = 0.972, r = 0.088 for n = 4. 

In the intermediate values of a between (0, oo) we need 
to solve Eq. (gDJ) for x by using Eq. ([21]). When n = 2 
the field value x can be expressed as Eq. (|23|) . in which 
case n s and r are known from Eqs. (|42[) and (|43|) for given 
values of a = m 2 /M 2 and iV. In Fig. [3] we plot the the- 
oretical values of n s and r for n = 2 as a function of cv 



(between 10 -8 < a < 10 8 ) with three different values of 
N (= 50,60,70), together with the la and 2a observa- 
tional contours constrained by the joint data analysis of 
WMAP7 @, BAO H3, and HST For a = lO" 8 

these observables are close to the values estimated by 
Eqs. (|44[) and (|45j) . whereas for cv = 10 8 they are close 
to those given by Eqs. (|46|) and (|4T|) . For larger a, r 
gets smaller whereas n s increases, so that the quadratic 
potential shows better compatibility with observations in 
the presence of the field derivative coupling to gravity. 

In Fig. [3] the theoretical values of n s and r are plot- 
ted for n = A as a function of a = X(M p \/M) 2 between 
(10~ 8 , 10 8 ) with N = 50, 60, 70. In the limit that a -> 
the quartic potential is outside the 2a observational con- 
tour for N smaller than 70. In the presence of the field 
derivative coupling to gravity the model can be compat- 
ible with the current observations due to the suppressed 
tcnsor-to-scalar ratio and the larger spectral index. For 

= 60 the parameter a is constrained to be 



a > 3 x 10~ 5 (95 % CL) . 



(48) 



and a > 4 x 10 -4 (68 % CL). For N = 50 the constraints 
are a > 2 x 10~ 4 (95 % CL) and a > 1 x 10~ 2 (68 % 
CL). 

Using the scalar power spectrum (f3"Tj) , the CMB nor- 



7 



malization by WMAP corresponds to 

2 



a 



M 



12^4 I Jf: I x eo 2 ( n + ax eo) - 2 - 4 x 10 ~ 9 ■ (49) 



where xqo is the value of x at N = 60. In the regime 
a > 1 one has x n+2 = [2N(n + 2) + n]n 2 /(2a), so that 
Eq. (gSJ) gi vcs 

— ~1.5xl0- 10 ^ for n = 2, (50) 
M p i M v ; 

-32 / -^pl 



A ~ 5.9xlO- 32 (^^J for n = 4. (51) 

In the case of the quartic potential it is possible to realize 
A ~ 0.1 for M ~ 2.8 x 10~ 8 M ph 



B. Exponential potentials 

Let us proceed to the exponential potential 

V{<j>) = V eP+/ M >* , (52) 

where Vo and /3 are constants. In this case one has 
e v = (3 2 /2 and r) v = (3 2 . Hence in standard slow-roll 
inflation we require the condition 1 <C 1. Moreover 
this corresponds to the power-law inflation without the 
graceful exit [23] ■ In the presence of the field deriva- 
tive coupling to gravity, however, the slow-roll parameter 
e = e v / A can be smaller than 1 even for steep exponen- 
tial potentials with (P > 1. Inflation ends when e grows 
to the order of unity, which is followed by reheating with 
gravitational particle production. This situation is analo- 
gous to that in brancworld inflation where the dominance 
of the density squared term (p ) can lead to cosmic ac- 
celeration for steep exponential potentials (42| . 

We focus on the case in which the condition 
V/(M 2 M 2 l ) 3> 1 is satisfied during the whole stage of 
inflation. Using Eq. ([T7| the field value <j>f at the end of 
inflation can be estimated as 



7 oe W//Mpi 
M 2 A/2j 



2 



(53) 



which implies that (3 2 3> 1. From Eq. (fl8|) it follows that 



y oe /30/A/ pl 



»4 



(54) 



In the regime ^4^1 the scalar spectral index (|3"2"j) and 
the tensor-to-scalar ratio (|3^|) reduce to 



(55) 
(56) 



2N + 1 ' 
16 



2N + 1 



which correspond to taking the limit n — > oo in Eqs. (|46[) 
and (|47)) . 



(a) N = 50 

(b) N = 60 

(c) N = 70 



0.10 



0.050 



0.0 




0.97 

n 5 



Figure 4: Observational constraints on inflation with the ex- 
ponential potential V{4>) = V r e w/M p 1 in the regime A » 1 
for JV = 50, 60, 70. The 1<t and 2a observational contours are 
the same as those in Fig. [2] Even the steep exponential po- 
tentials with [3 larger than 1 can be allowed from the current 
observational data. 



In Fig. H] we plot the theoretical values of n s and r for 
TV = 50,60,70 as well as the la and 2a observational 
contours. This shows that even the steep exponential 
potentials with j3 2 3> 1 are compatible with the current 
observational data. 

Using the WMAP normalization V s ~ 2.4 x 10" 9 at 
N = 60 in the regime V/{M 2 M 2 X ) > 1, we obtain the 
constraint 



/3~8.8x 10- 6 ^ 
M 



(57) 



For M/M p i < 10~ 5 one has j3 > 1. 



Natural inflation 



Natural inflation [24| is described by the potential 



V(cj>) = A 4 



1 + cos 



(58) 



where A and / are constants having the dimension of 
mass. In the absence of the field derivative coupling to 
gravity the above potential can be c omp atible with ob- 
servational data only for / > 3.5M p i [25|. Then a global 
symmetry associated with the pseudo-Nambu-Goldstone- 
boson is broken above the quantum gravity scale, in 
which case standard quantum field theory may not be ap- 
plicable. If the potential (|58[) originates from the string 
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axion, the regime / > M p \ is not generally realized [43| . 
This problem can be circumvented by ta king into account 
the field derivative coupling to gravity 1 [17| . 

In the following we focus on the case in which the con- 
dition A 3> 1 is satisfied during the whole stage of infla- 
tion (0 < 4> < 7r/). From Eq. (fl"7"| the end of inflation is 
characterized by 



14 



cosxf — — 1 
where Xf = 0/// an d 



VI67+I - 1 

47 



/ 2 A 4 

The number of e-foldings is given by 

N ~ 7 [^(X/) - -F(X)] , 
where % = 0// and 



(59) 



(60) 



(61) 



7-( X ) = 21n 



1 



1 



sin x tan x 



+ 2 In (sin x) + cos x • (62) 



The scalar spectral index and the tensor-to-scalar ratio 
are 



2 2 — cos y 
1 - - 



7 (1 + cosx) 2 ' 
1 — cos x 



7 (1 + cosx) 2 



(63) 
(64) 



For given 7 the field value Xf is known from Eq. (|59|) . 
Since x is determined by Eq. (|61|) for given values of 7 
and iV, we can numerically evaluate n s and r as functions 
of 7 for several different numbers of e-foldings. In Fig. [5] 
we plot n s and r in the range 7 < 7 < 10 6 , together with 
the lcr and 2cr observational contours. For increasing 7, 
both n s and r get larger. 

In the limit that 7 ^ 1 it is possible to estimate n s and 
r analytically. From Eq. (|59|) one has cosx/ — — l+l/^/y 
and hence Xf 1S close to tt, as (tt — X/) 4 — 4/7- Using 
Eq. dSU, it follows that iV~ (7/16)(tt-x) 4 -1/4. Then 
Eqs. (J63J) and (JM1) reduce to 



1 



AN+ 1 ' 
16 



47V + 1 



(65) 
(66) 



These values correspond to n = 2 in Eqs. (|4"6"]) and (f4"T|). 
This comes from the fact that in the regime 7 3> 1 in- 
flation occurs around the potential minimum at <f> = irf. 
As in the case of the quadratic potential V(<j>) = m 2 (f> 2 /2 



1 See Refs. [44j ] for other attempts to realize / < M p \ in large-field 
axion models. 



(a) N = 50 




(b) N = 60 




(c) N = 70 








//""< a) (b) 








































* * Larger 7 

1 
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Figure 5: Observational constraints on natural inflation with 
the potential V(cj>) — A 4 [1 + cos (</>//)] in the range 7 < 7 = 
f 2 A 4 /(M 2 Mp l ) < 10 6 with .4 > 1 for three different values 
of N. The la and 2a observational contours are the same as 
those in Fig. [2] For larger 7, both n s and r increase toward 
the values ra s = 1 - 6/(4iV + 1) and r = 16/(4AT + 1). 



with a = m 2 /M 2 3> 1, natural inflation with 7 ^> 1 is 
compatible with the current observations. 

For N — 60 the parameter 7 is constrained to be 



7 > 7.5 (95 % CL) . 



(67) 



and 7 > 9.8 (68 % CL). For N = 50 the constraints are 
7 > 7.8 (95 % CL) and 7 > 10.6 (68 % CL). In the regime 
7 3> 1 the WMAP normalization at TV = 60 gives 

/~ 8 x 10 4 7 1/4 M. (68) 

If M < 10- 5 7-V 4 M p i, then one has / < M pl . 

D. Hybrid inflation 

Finally we study hybrid inflation with the potential 



1 



V(4>) =V + ^-m 2 <f> 2 



(69) 



where Vq and m are constants. Inflation ends at a bifur- 
cation point given by <j> = <j) c due to the appearance of 
a tachyonic instability driven by another field X- As m 
the case of the original hybrid inflation [35| we focus on 
the regime Vq 3> m 2 <f> 2 /2. Note that in another regime 
Vq <C m 2 (f) 2 /2 the situation is similar to that in chaotic 
inflation discussed in Sec. IIV Al 
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Using Eq. (|18|) the field value can be estimated as 



: exp 



1 



-iV 



where fi and v are positive constants defined by 
Vo 



m 2 M 2 



M 2 Mpj 



From Eqs. (J32J) and ([39]) we obtain 



2v 



1 + fX 



8v 2 ( <j) c 

i + /i V^pi 



3 K-i)jv 



(70) 

(71) 

(72) 
(73) 



which mean that the scalar power spectrum is blued- 
tilted (n s > 1). Compared to the standard hybrid in- 
flation, the presence of the field derivative coupling to 
gravity (fi > 0) leads to n s close to 1 as well as the sup- 
pressed tensor-to-scalar ratio. In the limit that fi ^ 1 
one has n s — > 1 and r — > 0, i.e. the Harrison-Zcl'dovich 
(HZ) spectrum. The HZ spectrum is under the observa- 
tional pressure [f| , but this property is subject to change 
depen ding on the assumptions about the reionization sce- 



nario 



45j . The future high-precision observations will 



provide a more concrete answer about this issue. 

In the regime fi 1 and v <C 1 one has <j) ~ 4> c for 
N ~ 60. From the WMAP normalization it follows that 



^ 5 .3xl0- 4 

V (p c 



(74) 



If fx = 10 2 , z/ = 10 2 , C = 0.2M p i, for example, 
M ~ 10~ 8 GeV. In this case Eqs. (|72J} and ((731) give 
n s ~ 1.0002 and r ~ 3 x 10~ 7 , which is close to the HZ 
spectrum. 



V. CONCLUSIONS 

In this paper we have studied observational con- 
straints on a number of representative inflationary mod- 
els with a field derivative coupling to the Einstein ten- 
sor, i.e. G fJ, ^d fl cj)d v (j)/(2M 2 ). Such a non-minimal deriva- 
tive coupling has an asymptotic local shift symmetry 
for a slow-rolling scalar field satisfying the condition 
e = {d(j)) 2 /{M 2 M 2 X ) < 1. Since the strong coupling 
scale of the theory is around the Planck scale for e <C 1, 
quantum corrections to the inflaton potential can be sup- 
pressed during inflation. 

The non-minimal derivative coupling to gravity leads 
to a gravitationally enhanced friction for the scalar field. 
This property allows us to accommodate steep potentials 
with theoretically natural model parameters for realizing 
inflation. Not only the quartic potential V(<f>) = A<^ 4 /4 
with A ~ 0.1 but the potential V{<j>) = A 4 [l + cos(0//)] 



with / <C M p \ gives rise to cosmic acceleration con- 
sistent with the amplitude of the CMB temperature 

anisotropies. Moreover the exponential potential V{4>) = 
y oe /90/M pl; which often 

appears after the compactifica- 
tion of extra dimensions in higher-dimensional theories, 
can lead to inflation even for f3 2 ^> 1 in the presence of 
the field derivative coupling to gravity. 

For the potential V((f>) = (\/n)<fr n of chaotic in- 
flation the tensor-to-scalar ratio r decreases for larger 
a = AMpj~ 2 /Af 2 , whereas the scalar spectral index n s 
increases. In the limit that a — > oo, n s and r approach 
the values given in Eqs. (HHJ) and (|4T]) . As we see in 
Figs. [5] and [31 for both n = 2 and n = 4, the asymptotic 
values of n s and r are within the la observational bound 
derived by the joint data analysis of WMAP7, BAO, and 
HST. For the quartic potential with N = 60 we found 
that the parameter a = X(M p \/M) 2 is constrained to be 
a > 3 x 10~ 5 at the 95 % confidence level. 

For the exponential potential V{4>) = Vbe^/ Mpl 
the asymptotic values of n s and r in the regime 
U/(M 2 M 2 j) > 1 are given by Eqs. (J5SJ and ([56]>. with 
fi constrained as Eq. (|BT|) from the WMAP normaliza- 
tion. Figure 0] shows that steep exponential potentials 
with j3 2 ^> 1 can be compatible with the current obser- 
vational data. 

In natural inflation with the potential V(<f>) = 
A 4 [l + cos(0//)] the obscrvables can be parametrized 
by the parameter 7 = f 2 K A /(M 2 M^) in the regime 
U/(M 2 A/ 2 j) > 1. For larger 7, both n s and r tend to 
increase toward the values given in Eqs. (|55|) and (|66p . 
These asymptotic values correspond to those derived in 
Eqs. (gBJ) and (gTJ) for n = 2 with a > 1. This prop- 
erty comes from the fact that for 7 ^> 1 inflation occurs 
around the potential minimum at <f> = irf. The observa- 
tional bound on the parameter 7 is found to be 7 > 7.5 
for N = 60 (95 % CL). From the WMAP normaliza- 
tion the symmetry breaking scale / is constrained to be 
/ ~ 8 x 10 4 7 4 / 4 Af , which can be smaller than M p i for 
M < 10- 5 7- 1 / 4 A/ p i. 

In hybrid inflation with the potential V((f>) = Vq + 
m 2 <p 2 /2 (where Vq ^> m 2 (f> 2 /2) the field derivative cou- 
pling to gravity leads to the blue-tilted scalar power spec- 
trum close to n s = 1. Compared to standard hybrid in- 
flation, the power spectrum approaches the HZ one, i.e. 
n s = 1 and r = 0. The HZ spectrum is in tension with 
observations, but we have to caution that this property is 
affected by the assumption of the reionization scenario. 

If future observations such as PLANCK [46| can con- 
strain the tensor-to-scalar ratio at the level of r < 
O(0.01), this will allow us to place tighter bounds on 
the inflationary models with the field derivative coupling 
to gravity. We hope that we can discriminate between a 
host of inflationary models within the next few years. 
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